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For an infinitely long plate strip of constant curv- 
ature there is set up the buckling condition for combined 
shear and axial stresses. Assumptions are made of the 
vanishing of the moments acting on the longitudinal edges 
and of the vanishing of the tangential displacements u 
and the transverse stresses cy (longitudinal stif feners 
nonresistant to bending in the circumferential direction) 
The "buckling condition is evaluated for the case of pure 
shear and an approximate formula is derived for the criti- 
cal stresses for large curvatures. In connection with 
the discussion of the methods and results there is de- 
fined the range of validity of the "buckling condition set 
up. 



I, INTRODUCTION AND STATEMENT OF THE PROBLEM 



The stressed skin covering of fuselages and wings 
is sectioned off by longitudinal stif feners into panels 
of small curvature by these panels constituting curved 
plate strips. If a membrane type stress condition is 
produced in these structural parts the equilibrium con- 
figuration reaches the limit of stability at a definite 
load. Since each plate strip represents only a small 
portion of the stressed-skin structure the state of stress 
may be assumed as homogeneous before buckling sets in and 
there then arises the important problem of determining 
the critical loading limit for shear or compressive 
stress or for the two combined. 

The stability problems of a curved plate strip may be 



* "Die Stabili t&tsgrenze eines gekrtlmmten Plat tenstreif ens 
bei Beanspruchung durch Schub- und L&ngskr&f t e . " Luft- 
f ahrtf orschung, vol. 15, no. 10/ll f Oct, 10, 1938, pp, 
517-28. 



3 A. C .A . Technical Memorandum Ho. 898 



classified both as to the type of loading and the end 
* conditions. Since the resistance to buckling of the lon- 
gitudinal stif f ener s is considerably higher than that of 
the plate, a rigid support in the radial direction must 
be assumed at the edges (with the limiting cases of sim- 
ple (no moment) or clamped support). The tangential dis- 
placements, on account of the large extensional stiffness 
of the stiffeners as well as the radial displacements are 
to be set equal to zero at the edges. For the tangent to 
the plate at right angles to the edge, the following lim- 
iting conditions are possible: The vanishing of the nor- 
mal displacements (flexurally rigid edge stiffeners) or 
the vanishing of the normal stresses (edge stiffeners 
without f lexural resistance in the plate surface). 

A number of the problems here discussed have 
previously been treated. The critical stress in the case 
of compressive stresses has been determined by S. Timo- 
shenko (reference l) for an infinitely long plate strip 
with constant curvature under the assumption of simply 
supported longitudinal edges, Leggett (reference 2) has 
investigated the critical shear stress of a plate strip 
of small curvature for simply supported and clamped (flex- 
urally stiff) edges. His numerical evaluation of the 
buckling condition set up is valid, however, for only 
plates with very small curvature. In view of airplane 
construction requirements, however, his computations must 
be extended to larger curvatures . 

In the present investigation the buckling condition 
is set up for an infinitely long plate strip of constant 
curvature for the case of moment-free support at edges 
f or combined shear and axial loading . The buckling con- 
dition is evaluated for the. case of pure shear stress and 
for larger curvatures an approximate formula is derived. 
The evaluation of the buckling condition for the general 
case of combined shear and axial stresses will shortly be 
published. In contrast to what was done by Leggett we 
shall start with the complete fundamental equations for 
the displacements. With the aid of an expression by fiitz 
for the buckling functions the factors neglected by 
Leggett will then be justified for a certain range of 
plate dimensions . 
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-> I I ■ G-M32RAL PRINCIPLES 

Let there "be given an infinitely long plate strip of 
width 13, thickness s, and constant radius of curva- 
ture r (fig. 1). The zy;.y, z coordinates are taken re- 
spectively in the directions of the generators, along the 
circumference and along the radius, the x axis coinciding 
with one of the edges of the strip. Corresponding to the 
coordinates are the displacements u, v , w, of the middle 
surf ace, which carry the plate from the initial equilibrium 
condition to the neighboring condition on reaching the 
limit of stability. To produce these displacements an 
amount of work A^ in deforming the plate (potential of 
internal forces) and an amount of work A a against the 
external forces is required* Denoting the negative prod- 
uct of the external forces by the displacements of their 
points of application by V = -A a the change in the total 
potential energy due to the virtual displacements u, v, 
w, is equal to 

6 n a 4 + 7 (i) 

The linear components of the displacements u, v , w, in 
6 n must vanish since in the initial state the plate 
strip was in equilibrium. Below the stability limit the 
square components Q of 8 IT ar£ positive (minimum of 
potential energy in the stable equilibrium state), At 
the stability limit at least one set of u, v, w must be 
possible for which the square components Q vanish but 
none for which they are negative. From this minimum prop- 
erty it follows that for all admissible variations it is 
necessary that 

M=0 (2) 

(See reference 3.) 

To set up the expression for the deformation work it 
is necessary to find the relation between the plate de- 
formation and the values of the displacements. The work 
done in deforming the plate is made up of the work done 
by the basic stresses, uniformly distributed over the 
plate thickness, in deforming the plate middle surface 
and in the deformation work of the additional buckling 
stresses. In the first portion the deformation magnitudes 
enter linearly while in the second they enter squared. 
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Since the expression for the strain energy must also con- 
tain the square terms in the displacements u, v, w, the 
squares of the displacements must be taken into account 
in the deformations of the plate middle surface. 

Let t lf tgi 1 3 he the unit vectors in the direc- 
tion of the generators, the circumf erence, and the radius, 
respectively-, at any point of the middle surface, of. the 
plate (fig. 2). The displacement v will here he assumed 
along the tangent and not along the circumference. Let 
the position vector of the middle surface of the plate 
"before the deformation be R = S 0 . The square of the 
length of the line element before the deformation is then 



d s Q 2 




After the deformation 



S = B 



U ti + T tg + W tj 




(1 + u x ) tj + v x t 2 + w x t 3 



9 £ 

IT 



= R 



far - 



+ V 



y + J ) ** + ("y 




and 



ds a = dx s [l 



•f 2 u x + u z 2 4- r x 9 ' + w : 



x a 3 



+ 2 dx dy u. 



+ dy d j 1 + 2 



L 




For the changes in the coefficients of the square of the 
linear element of the middle surface of the plate we then 
have 
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a 



2 



? 22 = 2 (v + r ) + u r + K + ?)* * (*J " fY 

^ia = u y + v x + u x u y + v x ( v y + j ) + *x (*y - J ) 

The geometrical significance of these deformation magni- 
tudes is obtained from their relations with the exten- 
sional and angular strains of the surface. fe have (ref- 
erence 4) 



€ x = 7 

7 y = + %l - 1 = v y + I + § n y 3 + |(w y - l) 



sin 7 xy 



Jl + 7 X1 7l + 7 2 a + 7 *> + V> 



In the linear portions the magnitude 7 1S is equal to 



7 Y v ; values 



and 7 j 



differ from the c orre- 



ry, VJ .x^ v«,xu.«o J^S2 

sponding strains € x and c y by the factor 2. 

Por the deformations at the distance z from the 
middle surface there are obtained under the usual assump* 
tions of plate bending theory the following relations 



€y - Vy -h 



w 



r + 2 



- 2 



r -f z 



*yy 



r •+■ z 



'xy 



r + 



0 + — )• 

\ r + z / 



f aey 



(The detailed derivation of these and the following rela 
tions are found in reference 5.) Since these magnitudes 
enter only into the strain energy of the additional (buck 
ling) stresses only the linear portions need be taken in 
to account. By Hookers law the additional stresses in 
the plate after buckling are found to be 
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i - v 3 ' 



l _ v 



2(1 + u ) 



'xy 



(6) 



These equations, integrated over the plate thickness, tak- 
ing account of the trapezoidal form of the plate element 
in the cross section, give the following resultant forces 
and moments (fig. 3.) 



N„ = 



w xx 



D ( u x + uv y + v f )" | 

D (V-U X + Vy +.2) + | (Wyy + ) 



N 



xy 



1 - v , K 1 - v /v x 

— ^— <»y ■+ v x) + 7 —j- V^r 



2 

1 - V 



w xy ) 

/ \ K 1 - U / u v \ 



= K 



My = K 



M 



xy 



My x - K 



where D = 



(» 'XI + »yy + JT ) 



1 V 8 



is the axial rigidity of the plate and 



K « 



E s 3 



12(1 - u s ) 



the flexural rigidity. 



The deformation work, of the additional stresses per unit 
area of the middle surface ifc 
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| J ^xH + cr y cjr + t 7xy) (l + |) da 



+S/g 

a l = i J ((T x € x + dy€y + T 7 x y) (l + » ) dZ (8) 

-s/s 

Substituting equations (6) there is obtained 
+s /a 



a,-. I E 
2 1- us 



[(€ x +€y) S -2(l-u) (e x € y - lY xy )j (l+i) dz (8a) 



- S /2 

For the infinitely long plate atrip the displacements 
must be periodic along the longitudinal direction. We 
shall denote the half-wave length by I. for the strain 
energy of the additional stresses over the plate width b 
and a complete wave length 2 I there follows from (8a): 

2 I t> +s/s 

0 0 — S/ 2 

- 2(1 - v) (c x € y - -1 7 xy )j (l + f)dxdy dz (9) 

Ihe basic (original) stresses are made up of the membrane 
shear forces t s the longitudinal forces cr s, where 

positive values of cr denote tensile stresses/ Since the 
moment of the additional stresses about the normal to the 
plate 

r(N xy - S yx ) + Myx 



vanishes identically, as can be readily seen by substitut- 
ing the deformation values from equations (7) , the moments 
of the basic stresses at the deformed element must be also 
in equilibrium. (See reference 5, p. 200.) For the shear 
forces on the deformed shell element, we must therefore 
substitute ts(1 + F y ) dy and ts(1 + l x )dx. For the 

longitudinal forces, the assumption may be made that they 
act along the direction of the edges. The strain energy 
of the basic shear forces per unit area of surface (see 
fig. 4) is obtained as 
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a 2 = T s ( 1 + T x ) (1 + T y ) sin 7 xy (10) 

With the aid of the third of equations (4) there is ob- 
tained 

a 3 - rs 7 13 (10a) 

The strain energy of the longitudinal forces per unit 
area is obtained as 

a 3 = a s 7 X (11) 

Integrating over the plate width and a wave length, there 
is obtained from equations (10a) and (11) for the strain 
energy of the basic stresses 



a I b 

. a r 



A 3 + A 3 = TS 



y J [ u y + v x + u x% + ^x (vy + J) 



( w y - 7 )J dx dy * °" s J I Ux * ^ v x 2 + w x S ) j dx d y 



0 0 

2 t "fa 

n 

/ v \ i_ 

+ W X 

'° ~* (12) 
The potential of the external forces is 

si b 

' • / ' v ~- u J — U ' 

'0 



«5 * U 

= - T.s / (u y _ b - u y==0 )dx - ts I (v x=2l -v x==0 )dy 



b 



'Jo 



As may be easily seen from equations (12) and (13) the 
linear strain terms in the sum cancel each other (condi- 
tion* of equilibrium). In the change in potential energy 

6 n = Ai + 7 = Ai4 A 2 + A 3 + T 

there remain only the quadratic portions Q. By substi- 
tuting the displacements (5) in equation (9), there is 
obtained 
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a I b 
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0 0 a t Is 



0 0 



-3 (u x + Ty + - 2(1 - v) [(•»„ w yy - w iy s ) 

(Vy + Z )+ w x (w y - l)J dx dy 



alb 



+ T S 



in 



U X Vy + V X 



0 0 

2 I b 



+ as X J J (v x s * w x 2 ) dx dy 



(14) 



O 0 



for simply supported 'longitudinal edges w == 0 and 
Wyy = 0 at y == 0 and y - b. Furthermore, the tan- 
gential displacement u must also vanish at these edges. 
Making use of these end conditions and the periodicity 
of the displacements along the longitudinal direction 
and integrating in the first term of the sum partially 
with respect to x and the second term partially with re- 
spect to y, it is readily seen that the portions 

u x v y - u y v x> w xx*yy - .*xy* , and v x 23. - v y 



do not contribute to the strain energy. From the minimum 
condition §Q » 0 we then have 
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-5^ = 



0 o 



u xx ■+ g u yy + ~ 



g-^xy + v ~ 



+ 2 1 i u 



8 u dx dy = 0 



(15.,) 
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3 I 



} 



dx = 0 (15. J 



2 I D 
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) 



where 



£ 1/8 X 65 TS /. n<l2 x T 



w xx] 



(15. 3 ) 



and 
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Since the equations (15) thus set up must vanish f or 
each variation 5u t 6v, and dw, the expressions in 
hr.ackefrs under the integral signs are to "be set equal to 
zero. This leads to three linear differential equations 
for the displacements u, v t w f which equations may also 
he derived from a consideration of the stahility of a 
shell element. For longitudinal memh e r s that are flexur- 
ally yielding in the circumferential direction Sv is 
not equal to zero at y = 0 and y= h and from the 
"boundary integral of equation (I5 a ), there is then oh ~ 
tained as the boundary condition for the displacement at 
y = 0 and y = h . 

v y + q x v x = 0 

Considering feqixation (7. 3 ), . it may he seen that the 
normal stress and the component of the shear force in the 
circumferential direction at the deformed edges y = 0 
and y - h must he in equilibrium* . The complete "boundary 
conditions at y = 0 and y == b are therefore 



w = 



f ¥7 



v y + V* 



u = 



v_ = 



0 
0 
0 
0 



In an exact determination of the crit 
the solutions of the differential equation 
he substituted in the "boundary conditions 
then obtained eight homogeneous equations 
edge) for u, v, and w. The problem is • 
smallest value of q x for given q s (or 
which a nonidentically vanishing solution 
placements is possible. The differential 
he solved by assuming for the displacement 
sions 



ical stresses, 
s in (15) must 
(16). There are 
(four at each 
then to find the 
conversely) for 
f or • the dis- 
equations may 
s the expres~ 



u 



= A sin i (Xx - 



« B sin i ( 



Xx - 



W == 0 COS i (XX 

r 



m y) 
m y) 
m y) 
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The continuation of the method leads, however, to compu- 
tational difficulties which, like those met in the case 
of a twisted uns tiff ened circular cylindrical shell of f i- 
nite length (reference 5, p. 203), exceed reasonable 
bounds. The computation will here be carried out with the 
aid of a Hit is expression in the form already employed by 
Leggett for the strip of small curvature and by U. Seydel 
for the flat strip (reference 6). 



III. SITTING UP OF THE BUCKLING CONDITION 
WITH THE AID OF A RITZ SUBSTITUTION 



The Ritz method consists in approximating the dis- 
placements by using a combination of partial expressions 
in which undetermined values enter linearly, each expres- 
sion satisfying the boundary conditions. The values are 
then so determined that the approximating expression, sub- 
stituted in equation (14) , reduces the latter to a mini- 
mum. The calculus of variations problem then goes over 
into an ordinary minimum problem. If the minimum condi- 
tions for the quadratic terms in the change of potential 
energy are set up before the integration (method of 
Galerkin, reference 7) there is obtained again the system 
of equations (15) where, for example, there is to be sub- 
stituted 6 u = 2 ~— 6 A m where A m are the undeter- 
m 3 A m m m 

mined values entering in the expression for u. Since 
only known . functions now stand in the integrands, the in- 
tegration can be performed and on the left-hand sides 
of the equations, there are then obtained polynomials in 
the variations 5 A m , etc., the undetermined values them- 
selves occurring in the coefficients'. Since the equations 
must be satisfied for each system of values 6 A m , the 
coefficients must vanish separately and this leads to a 
system of equations for computing the undetermined values. 
This system of equations will have a solution for the 
displacements other than zero only if the determinant of 
its coefficient's vanishes. There is thus obtained the re- 
quired condition for the reaching of the stability limit . 

The oblique buckling waves that are obtained in the 
presence of shear are represented in the expressions for 
the displacements by two partial waves for each displace- 
ment shifted in phase by a half -wave length in the longi- 
tudinal direction: 
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U ~ COS 



TT X 



2 A m sin 



m tt y 



it x v a • m tt y 
- sin -~r- E A orn sin ~™.~- z L 
I m 2m ^ 



v =5 sin 



TT X 
I 



w == sm 



w ^im 



cos 



m tt y 



■+ COB 



TT x 



m 



B 2m cos 



m tt y 



tt x n . in tt y 

-r— 2 Cim sin — ~- 



m 



+ C os H-S S Cam sin SJLZ 



m 



V — -(17. 



with m == 1, 2, 3 . . . • The boundary conditions (16. 1 _ 3 ) 
are exactly satisfied by the chosen expressions, while 
"boundary condition (16. 4 ) is satisfied except for the term 
q v x , As may be seen from equation (17. a ) the derivatives 
ox v with respect to x and y for the range (0 < y < 
b) are of the same order of magnitude. Since q x = 

X (1 - v B ) <^ 1, it follows that an exact satisfying of 
E 

the boundary condition (16. 4 ) would require in the expres- 
sion chosen for v only a slight correction which is dis- 
pensed with here. 



As will be shown below, in the first partial waves we 
require only the even terms in m and in the second only 
the odd terms (or conversely). The number of terms taken 
gives the degree of the approximation for the buckling 
functions and for the critical stresses. As the computa- 
tion for the case of pure shear shows, the series of ap- 
proximations for the critical stresses converges very rap- 
idly. The summation for the partial waves in the trans- 
verse direction of the plate strip needs to be extended 
only to m ,== 4. 

The functions thus chosen for the displacements u, 
v, w are now substituted in equations (15 vl _ 3 ). Since 
the boundary condition (equation 16. 4 ) is not exactly sat- 
isfied, the boundary integral in equation (15. 2 ) will 
yield a Small contribution. Making use of 
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(18) 



there is obtained the following 
system of equations s 

+ »!«/»- C,»^»==0 
«» t- ft 0 2 m ^r^i 

: + Bznfn — Czn — gn^O 
Aw fn + BinK + 9x P 2/ 5am T n 8 ~- ^ a 



5 im " 



m 

* w ji r • 1 n tit m n* — to* 
tit 

.8 6 



1 tt 8 4-TO» 
» 8 — TO 8 

TO 



b D & * fl 8 6 x? 1 to » 

- A t „£-g„~ B,„^ fc„ + ft ^f^X *i« jpz^ 



where, for briefness, there has 
been set 



(22) 
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In the summation terms m is to he taken so that 
n * m is an odd number (see equation (18)); i.e., for n 
even, m must he odd and conversely. The structure of 
the system of equations is shown by the scheme of figure 
5 which represents the part of equations (21, -^)lritb the 



unknowns C ln , 0 sn .. It may he seen that the entire sys- 
tem breaks up into two independent partial systems whose 
coefficients are distinguished hy the heavy and light 
dashes, from equations (21) it is further evident that 
the coefficients of 0 ln in 6 qua t ion (21. x ) are equal to 
those of C gn in equation (21. a ) whereas the right upper 
region in figure 5 and,, the lower left region contain equal 
and opposite coefficients. If therefore the unknowns 
0 V1 » 0 13 . . • , C 2S , C g4 *.. are computed from one partial 

system (light dashes) the solution of the other partial 
system can he found from the relations 



12 



'32 



'21 



C !4 ~ G 24 



°33 ~ -°13 



- -- (23) 



The same holds true for the A's and B's since the en- 
tire system (19-21) is constructed in exactly the same 
manner as the portion shown in figure 5. 

Substituting the equations (23) and the correspond- 
ing expressions f or A and B in (1?) , it . is seen that 
the two partial systems of figure 5 represent the same 
buckling form except for a phase displacement of 1/2 so 
that it is sufficient to investigate one of them in which 
the sum of the two subscripts of the' unknown is odd. The 
pairs of equations (19) (20) and (21) can then each be 
combined into one as follows: * ". 
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A n a n 



n m 



n 3 - m 2 



t-D m 



+ B n f n - 



°n rr six .= 0 



x tt m 2 n 2 »- m 2 



m 



b 8 " b ■ 

°n rr h n + a, P - ~r g o m . ~ _ (-i) m = o 



13 b 

L n — Sa - 3n -~- h n + 

TFT TTT 

8 t> 



n 



+ qp-~2 B m — _ 

i tt trr m n 2 - m 2 



(-l) m 



+ C 



n 



rr 3 r 2 



c n - q 6 _ 2 C m 



n m 



i tt m n 2 m 2 



(~D m - 0 



(24.^ 



where n and. 



m 



now run through both even, and odd values 



For n = 1, 2 ... i and m - 1, 2 ... i (in each 
t«erm n ± m is odd) , there is obtained from the coeff i- 
cients of the system- of equations (24 , x „ 3 ) a determinant 
symmetric with respect ot its main diagonal, of . the 

i^d. order. The required buckling condition is found in 
the following manner* Let us consider the first two sys- 
tems of equations of (24)'. for a finite number of terms 
n = 1, 2. • . • i the undetermined* values A and B can 
be expressed in terms of the' C 1 s by Cramer 1 s rule. From 
the consideration of the values a n , b n , f n , g n , h n , and 

the summation terms, it may be seen that the determinant 
of the denominator is of the order of magnitude i ! 4 . and 
similarly the determinant of the numerator of A n 



num er a tor d e t er m i nan t of B 

of magnitude i i ! 4 , 
n 



The 

however, is of the order 



By substituting the obtained expres- 
sions in the system of equation (24. 3 ) there is obtained 
a determinant of the coefficients of the i th order. By 
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dividing through by the term k IIeI + 0 s ) 

b 4 

t er m i nant may "be b r o ugh t into a 11 or ma 1 form; 



this' de- 



1 + c 



22 



11 13 
1 + C 

0 32 
• • * • • • 

where the infinite series 



21 



3 1 



13 
0 33 

1 * C 33 
• • • *.« 



• » ♦ 



• • * 



* • • 



(25) 



C ll + C 12 + C 21 + G 13 + C 22 + C 3i + 



(26) 



converges absolutely, since n~ 2 is the highest power of 
n occurring in c a p- This, however, is a suffi- 
cient condition for the convergence of the determinant 
(25). The latter set equal to sero gives the buckling 
condition for the curved plate atrip under combined shear 
and axial stresses* 

In the elimination of A and B it is necessary to 
solve determinants of the 12 order. With certain re- 
strictions for the plate dimensions, which, however are 
within the requirements of airplane construction, the com.* 
putatipns may be considerably simplified; Equations 
(24, 1 _ 3 ) are first transformed as follows: 



A n f a n - U i n3 n m ma (-D m V B tt ffc- C n JL. g n 
A n f n + B n U n - x 5 |a I * s t * 8 <-l) m 1 

nr I . m 0 n n a _ a a J 



= 0 



= 0 



(-1)™ 



* - \ 2 £a 

-n - a m 7T 



n m 



•&n n 2 ~ m s 



-J. 



v <27. ) 
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where, for "briefness, there has teen set X = <li 3 ■■■«%= In- 

TT 

troducing the values 

4 tt 2 K tt2 I s 2 



p* = 



3 1 - u 2 b s 



b a s 

as reference values for the critical stresses and remem- 
bering that ■ q x s= -L (1 - tf 2 ), there is obtained for X 
the expression ... 

T ft tt e 2 

(38) 



X = - p — 1 



2 



The reference value p* thus introduced is the critical 
compressive stress of a flat plate of width b for sim- 
ply supported longitudinal edges. For the flat plate 
t/p* ~ 1.334 and £ == 0.793; for a curved plate t /p* 

increases with the curvature but the product X- £ re- 

p* 

mains, as further computation will show, approximately 
constant* Thus X on account of s/b <c 1 is also a 
very small magnitude. 

For the f lat plate (r = co) equations (27 . lmm& ) sepa- 
rate from equation (27. 3 ) .. It is known that the vanishing 
of the brackets after G n in equation (27. 3 ) gives the 
smallest value for X and that An ^ Bn 5 0. Also in the 
case of plates with small curvature it is assumed in a 
first approximation that the X term in the bracket after 
On equation (27. ') remains the essential buckling term of 
the system of equations. If the X terms in equations 
(27 . 1 _ 2 ) and the small magnitudes k and q 3 are dropped 
then A n and B n may be expressed in terms of C n from 
these two equations, as follows: 



A n = 



(n 3 + p»)« 

n g + (3 + v) ff a 
(n 3 + p 3 ) 2 



(29) 



from which is obtained: 
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^2L 
An 



m 3 - v p s (a 3 t g a ) s O m 
n 2 - "u j3 8 (m a + p a ) a °n 



B m m m 3 + (2 + v) 08 («• + ps>a o m 
B a n n 3 + (2 + v) pa ( m s + ft 8 ) 2 C n 



(30) 



With the aid of these first approximations the ratios be- 
tween the A terms and the parts a n , "b n , etc., in the 
coefficients will be set up as a check on the approxima- 
tions made. If jU 2 is neglected compared with n s or 
m2 and in addition there are neglected (for r 4 s °°) 
small magnitudes k and q. a , the required ratios are 
found to he the following: 



for A n in (27. x ) 



1 - V 



m 



For B 



n 



m O n n 2 - m 2 m 
in (27. a ) 
C 



£ (-D 



m 



m O n n s _ m s 2 Am n' 



For 



in (27..) 



G m 



X £ 3 



m 



For B 



n 



m O n n a - m a n 
in (27. 3 ) 



- (-l) a : 1 



»' (31.,-s) 



n 



1 

m 0 „ n s - • m s m 



1 (-1)° :. 1 



For 



in (27.,) 



m 



m °n n a _ m a 



n m 



T T 3 r 3 
b s 



(-1) 



?or r — > cd ratio (31. 5 ) will become infinite ♦ As shown 
by equations (32), however, the Jc and q s term in c n 
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is not to "be neglected compared to 1, but conversely, and 

there is obtained in place of (31. 5) ' 

* a r a" tt a 3 3 1 fl • (-D 13 : 1 (32) 

12 # 3 b 2 * p* H 

where the reference value p* has also been introduced 
in the term with , q g . * The ratios (31 * ^4) are independent 

of the width of the plate. For sufficiently small curva- 

turo ^''' ^g » l^j these ratios are smaller than {31. 5 ) by 

the factor tT 2 r2/b s and in addition smaller than (32) by 
n 2 or m 2 in passing to the flat plate on account of 

s 3 

— — <C 1. There has thus been confirmed the assumption 
b 2 

made above that the X term af ter G n in equation (27. 3 ) 

remains the essential buckling term also for the case of 
plate strips of small curvature. The magnitudes of the 
ratios thus set up are dependent of G a / C n . For values 

of n to which correspond large values of G n (in the 
present case for n ■=* 1 to 3, as will be shown with the 
aid of examples); i.e., for considerable portions of the 
system of equations the ratios ' (31. !- 4 ) are small and the 
first approximations for A n and B n , according to equa*- 
t ions (29) are sufficiently accurate. For larger values 
of n, that is, smaller values of C n , the accuracy 
will decrease but the corresponding buckling terms,. as 
will be shown on evaluating the buckling determinant, are 
of small significance. 

The determinant of the coefficients of (27) must be 
equal to zero for nonvanishing values pf A n , B n , and C nt 

If, in the course of the computation, the powers and products 
of k, q 2 and of terms with X are neglected and for 
A m /A n and Bja/B n their first approximations from (30) 
are substituted, we arrive at the following system of 
equations: 
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0 n {(1 - V a ) 3 4 + k f(n 2 + P 3 ) 4 

- 2 -iL-(n 6 + (4 - v) n 4 B a + 3 n 3 8 4 + P 8 s ) .... 

+ -4^7 (n 4 + 2(2 - V) n a 8 a + (4 - 3 V s ) 8 4 ) | 
tt r J 

.+ q a ?2£ 8* [(n a + 8 a ) a + ~ (n a + 2(1 + P) 8 a )]}' 

. qi ^^ (n8 + , 3 ) Sm2Cm ^(-ir {l 

b a 1 f m a - V 8 s a 

+ — 1 -- 6 

TT a r2 (m 3 + 8 2 ) 2 in 8 - v 8 2 



m 3 + (2 + v) 8 2 _ . » o2N m s + n 2 

+ ■ ■-■ ■ ■ . ' (n 2 + 2(1 + v) 8 s ) ■■■.„.„..,,■ 

n s + (2 + p) §2 2 n a 

- (m3 * P 3)S (n a + (2 + v) 6 3 > 
(n 2 +, p a > 3 

- S * + <g + (n a + (2 + v) 6 3 ) 1} ■ (33) 
n s + (2 + v) 8 s J 

The convergence of the determinant of the coefficient's in 
equations (33) may he easily proved if the equations are 

divided through hy k 17 (n 3 + j3 s ) 4 . There is then oh- 

h 4 

tained the normal form (25) with* the converging series . 
(26). The left side of equation (33) is exact except for 
the neglected products and powers of k and q s and has 
the same form as for a circular, cylindrical shell under 
compression and without stiffeners.* 

* (See reference 5, p. 196.) The circular, cylindrical 
shell under compression during "buckling is divided up by 
the nodal lines of the "buckled form in the x direct icin into 
curved plate strips with w ~ Wyy == u = Vy ~ 0 at the 
nodal lines; i.e., the same "boundary conditions 1 * which are 
considered here. The exact computations of W. Jlugge can 
therefore he used for the curved plate strip. 
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If we neglected the shear load- terms in equations (27) 
with the exception of the important "buckling term for small 
curvature strips at 0 n in the third equation, there would 
remain on the right side of equations (33) only the 1 In 
the "braces. Let us denote the other terms as correction 
terms and expand them, since £ a is a small magnitude, in 
powers of p. We "break off at the quadratic terms in £ 
and obtain for the braces the following expression 

f n b 2 Pi /I 1 \ 1/1 1 \ a a 1 *T> , aS 

V - LI fa + ^ + I fa " W " P " 35? p Ji (34) 

The above elimination process is suff iciently 'accurate 
if the correction terms (34) are small as compared with 1. 
For a plate strip with a central angle of 45°, i.e., with 

= -k corresponding to a circular, cylindrical shell 
Trr 4 

with 8 longitudinal st if f ener s , the correction terms for 
n ■== 1 and m = 2 (n ± m odd) on neglecting the small 
terms with amount to about -4 percent. On the left 

side of equation (33) the coefficients of k and q 2 
(square brackets) are arranged in powers of b /rtr . For 

b 1 4 

^ =5 j the portions with b '/rr 4 r 4 at k are negligibly 

small and the squared portions in h/nr amount to about 
12 percent of (n 2 + P 2 ) 4 for n == 1. For larger values 
of n these terms on both sides of equation (33) become 
relatively still smaller. Considering that the extension- 
al part of the deformation (1 ~ U 2 ) |3 4 on the left side 
of equation (33) and the bending portion with 

k TT 4 r 4 for n = 1 are approximately equal, as will be 
b 4 

confirmed in the evaluation for the case of pure shear, 
there is found for the term with b 2 /fT 2 r 3 in the : coeff i^" - 
cient of k only about -6 percent of the no-load terms. 
Hence, for pure shear load (q 2 = 0) the terms with 
b s /TT 3 r s , since they are subtracted on bath sides of equa- 

b 1 

t ion (33), may be neglected even for ~- — . 

-nrr 4 

In the coefficients of the longitudinal stress q g 
on the left side of equation (33) the squared terms in 

b /tjt for b/rrr = i amount to about +6 percent of (n s '+ p s ) s 
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for n ~ X. For pur e compression the dropping of ilxe 
terms with b 2 /7T a r 2 and b 4 /n 4 r 4 in the coefficient s 
of k and q g f since they occur with different ki^ns, 
would lead to an error for the buckling stress of about 
+13 percent. In the table given below the err or $ which 
would arise from neglecting the above-mentioned terms, are 
given as a function of b/rrr f or k =s X fXZ /t% s= 10*? 
and v - 1/6 1 (the exact values are taken from W. Flugge) : 







1 


1 


1 


1 1 


1 


1 


Ttr 


10 


8 


6 


5 


4 


3 


2 


Error in 
















percent 


1.4 


3.2 


6 


8 


13 


25 


64 



Allowing an error of 6, percent, the simplified computation 
may be employed up to central angles of 30° (b/irr == 1 /6) 
for the case of pure compression. 



Under the restrictions of 45° for the central angle 
for pure shear and 30° for pur© compression, there is ob- 
tained from equation (33) after division by 

k (n 3 + j3 3 ) 3 ^ n £ the following system of equations: 



rr 



(0 n n) 



TT 



32 p 



(n 2 + j3 2 ) 3 + 



(n a + p 2 ) 



+ 1 4 p 8 

2 \ 3 -p* 



n s - m 2 



(-l) ra * 0 



(35) 



where 



b 4 r 5 

oj = 12 (1 - v 2 ) ~— ~ 

rr 4 r 4 s a 



(36) 



x The buckling formulas with the above neglected terms for 
the case of pure compression become the following. (See 
reference 1.) 



- JL = i + » 
p* 



.a _ 1 / — 
\p* ~ 2 



16 



'or uu $ 16 



for a) 2 16 
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The two parameters of a curved plate strip . b /rrv and b/s 
still remain only in 0 and in the . refer ence magnitude p* 
From this it may he seen that the critical stress in the 
case of pure shear or pure compression may each be repre- 
sented by a curve for a plate strip with small curvature 
and in the case of combined stress by a single parameter 
family of curves. 

If, for briefness, we set 



TT 



32 § L 



(n 2 + r ) S +U) 



,4 



cr 2 

, + - 4 p 3 

(n 2 + p 3 ) 3 p* 



— (37) 



n 



then for a given a/p* the critical shearing stress is 
given by setting equal to zero the following determinant: 



A = 



3 p* 



15 p* 



I -I 

3 p* 



A — — 0 -A. 



1 T 
5 p* 



21 p* 



1 _T 

5 p* 



1 -I 

7 p* 



1_ X 

15 p* 



- A -A. 0 



1 J_ 

7 p* 



I -I 

9 p* 



21 p* 



1 jr 
9 p* 



— A ■- ■ . i n T 



^0 



(38) 



The dropping of the terms with, b 2 /tt 3 t q and b 4 /TT 4 r 4 



in the brackets after k 



TT 4 r 4 
b 4 



i n e qua t i on (33) amount s 



to neglecting the trapezoidal shape of the plate element 
in cross section. From equations (5. ) and (7, X ^ Q ) 
there are then obtained the following simplified equations 
for the plate strains and for the section forces: 
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N 



€ y - T y + £ - z w 



yy 



^xy = % + *x - 2 2 w xy ^ 

D (P + Vy + .5 ) 



xy = H yx * D (u y 4- v x ) 

My - I(V W„ + Wyy) 

M xy = M yx * X<1 - ^)w xy 



, (40. 



The curvature of the plate element in these magnitudes is 
expressed only through the term w /r J ' Of . the load terms 
there remain after striking out the portions with 

the buckling term -2 q x r w xy and -q 2 r " w x3C in 
equation (15. , a ) and. -q x p» : 2 C m ^ (-l) m and 



TT a r s 



w a m a * xn 2 



q a £ a in equation (24. 3 ), i.e., only the buckling 

terms which would remain in passing, to the flat plate-. 
In the equilibrium condition there is added to the por- 
tions for the flat plate the radial components of the 
ring stresses . ?he system of differential equations*' : r * 
(15. 1-3 ) now assumes the following form: 

2 yy 2 ^y r 



1 ~ V 

v -4- a ■ nuk tr 

yy . 



XX 



+ 1 t v u xv + .2*. = 0 
. 2 xy r 



p ( u u x + v y + + k r s Aiw - 2 q l ^ x y 



- = 0 
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By introducing a stress function, equations- (4-0, 1-2 ) can 
"be further combined into one. The equilibrium condition 
of the forces in the tangent plane to the plate element 
are now simply 

3 N x 3 JT XV 
3x d y 

9 y 9 x 

These equations are exactly satisfied "by the expression 
for the stress function 



From equations (40) there is obtained by elimination of 
the displacements u and v the relation between the 
stress function 0 and the radial displacement w 

I A A $ = w xx (42.,) 

The differential equation of (41 . 3 ) assumes the form 

$ xx + k rS A A w - 2 q i w Iy - q 2 w x2C = 0 (42. g ) 

The system of differential equations (41 #1 _ 3 ) or (42, x _ 2 ) , 
the validity of which was proved here for a central angle 
of 45° for the case of pure shear and about 30° for the 
case of pure compression, is the starting point in the in- 
vestigation of Leggett. 

17. EVALUATION OF THE BUCKLING CONDITION AND SETTING-UP 
OF AN APPROXIMATION FORMULA FOR THE CASE OF PURE SHEAR 



In the diagonal terms of the buckling determinant 
(38), the longitudinal stress <j x is to be set equal to 
to zero. Using in the expression values for n up to 
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n * 2, 3 etc., there are obtained from equations (38) 

for the determination of the critical shear stress the 
following approximation formulas: 

First approximation: 

(~r) a = 9 * % C43) 



Second approximation; 

T„ 



(yf = -^i^ (44) 

5 3 X l + 3 2 A 3 

Third approximation; 

V' ^3X7 5X15 / Vp*/ L 3 s 5 s 7 s 15 s J 

+ T x T 2 T 3 T 4 - 0 , ; (45) 



The further approximations need not be written out since 
the critical shear stresses hardly differ in the 
second and third approximation. In the diagonal terms T n 
there enters the ratio of the plate width to the plate 
wave length .0 ~ b/l as an unknown. The latter must he 
so chosen as to give the least critical stresses. 

In figures 6 and 7 the results of the computations 

are plotted against VJJ - Vl2(l - < ~- / Keeping 

rrr v s 

• the ratio b/s fixed, the curves show the dependence of 
the critical shear stress and the reciprocal wave length 
- -on the plate curvature. The first approximation for the 
critical shear stresses gives in the entire range values 
about 5 percent too high while the second and third prac- 
tically coincide. Por u) - 0(r == <») , there is already 
obtained from the third approximation the exactly computed 
value f/p* - 1.334 of Southwell and Skan (reference 8). 
figures 6 and 7 also contain the curves computed by Leggett 
for the other limiting case of a plate strip with flexur- 
ally rigid edge stiff eners (v ~ 0). The critical shear 
stresses lie considerably higher than in our ease of. a 
plate strip with longitudinal stiff eners flexurally yield- 
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ing in the circumferential direction == JO.) ... in prac- 

tical applications the edge condition will lie hetween 
these two limiting conditions and the assumption here mad© 
is on the safe side. 

From the consideration of equations (44) and (45) it 
may he seen that on dropping all terms except those with 
3 2 in the denominator there is again obtained the first 
approximation (43)* Since the latter differs from the 
further ones hy only 5 percent the other terms in the sec- 
ond, third, etc., approximations which are obtained using 
greater values of n in the assumed expressions should he 
very small, There is thus proved the statement made in 
section III that the errors made with increasing n in 
the elimination process are insignificant in the final 
result . 

As may he seen from figure 6, the critical-stress 
curve with increasing id approaches a straight line asymp- 
totically. The equation of this asymptote may he easily 
derived from the second approximation hy introducing a 
new unknown and striking out small terms. The individual 
terms in equation (44) are the following: 

T x « f(l + p 8 )* + a) 

32 P L 



T 2 = _2L- [(4 + p s ) 2 + 
32 0 ^ 




* I (9 + P 2 ) S + 



s 4 


(1 


+ > a ) 2 






(4 


+ p»)» 






(9 


+ |3 2 ) 3 



As a new variable there is taken the ratio of the h end ing 
portion (1 + p 2 ) S and the extensional portion a> ft ... 

in JO and set equal to z*. 



(1 + p*) a = z 



U) 



e 4 



After Borne computation, there is obtained 
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> 2 - * yST p + 1 - 0 



As is easily shown by the numerical computation, for a 
large value of u> the small root must be taken for (3, 
Neglecting small terms, there is obtained 



2 VuT 



TT 4 
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V U) 9 z, etc. 



(46) 



and substituting in equation (44) 



— * ^r) 



( 



35 + ^ 



a (t fa* ) 2 

From the minimum condition - — • ~ ■ = 0, ..there follows 

d(z s ) 



z = 0. 956 



1 . 045 

4 y— — 

V 0) 



(47) 



— ~ 0.30 Vo) 



p* 



(48) 



The product . — p remains for large value of u) constant 
P* 

and the extensional and- bending portions in $i on ac- 
count of z = 0.956 are almost equal. 



Prom figure 6 it may be ssen that the curve of crit- 
ical stress may be replaced by its asymptote for 
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V = 2.5, i.e., for <.u tt 40- with sufficient accuracy. 
The equation of the asymptote ' (48) for u>>40 may "be 
assumed as the approximation formula for the computation 
of the critical shear stress . After some computation 
there is obtained the .following "buckling formula 



T = 1.56 
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(1 ~ » 2 > 
hence with v = 0.3 



s /s 



3/4 "b v r 



T = 1.67 



s /IT 

W r 




(49) 



4 : . 

f or - / - < 4,3 or a) = 12(1 - ~ — - — 



s 



TT 4 r 4 s 2 



< 40, the 



critical shear stresses must be taken from the curve in 
f i gur e 6 , 



7, SXiMPLSS 



Ihe accuracy of the eliminati on process described in 
section III for a central angle of 45° 



( - = - ) 

\ -rrr 4 J 



which was assumed as the limit for the compu- 



tation carried out here, will be tested on some examples. 
There will first "be considered a plate strip with b/s = 
50. For u = 0.3 there is obtained from equation (.3.6) 
to = 172.5 and further f or j3 a = 0.1, p* = 2. 95. (See 
fig. \6.) From the buckling determinant (38) the coeffi- 
cients in the expression are computed to be 



0 3 = -0.0558 -G 



C 4 - -0.0065 O x 



G s = -0.413 Gj 
The buckling function then assumes the following form: 



— — w = sin' 



/0 - 1 B * [-0.413 sin UU. - 0.0065 sin -JLZ 
b l b b 



VOTT TT X f TT y „ „ . 3 TT 7 

+ cos - sin — - - 0.0556 sin - 

b L b, D • 
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In figure 8 the corresponding buckling surface is reprew 
sented means of contour lines. In figure 9 this is 
compared with the buckling surface of a flat plate taken 
from the paper by Southwell and Skan (reference 8). JPrpm 
a comparison of the two figures it may be seen that the 
wave length of the curved plate strip is larger and the 
folding angle made with the longitudinal axis smaller than 
for- the flat plate. 

The accuracy of the fir st approximations for A n and 
B n (equation (29)) depends on the ratios of the neglected 
portions in the coefficients of A n , B n , and C n , to the 
remaining portions a n , b n , etc*, in equations (31 t 1 _ 3 ) • 

For the absolute values of the ratios in equation (8I.«i) 
there is obtained 

for n » It 0.000615 
n ~ 2: 0.0145 
n = 3; 0.0198 
n = 4; 0.347 

The ratios in equations (31. 3 ) and (31. 3 ) are, as may also 
be seen from the form of the equations, somewhat greater 
for n = 1 than for the equation (31 1) but for n = 2 to 
4, smaller, from the computed values it may be seen that 
the first approximations for A n and *B n for n =1 to 
3 are sufficiently accurate. The larger error for n = 4 
is unimportant for the computation of the critical stresses 
on account of the small magnitude of the corresponding 
load terms. 



For larger values of u) (i.e., for larger b/s at 
the same central angle of 45°) the equations (46), (47), 
and (48) may be substituted for the values T n , (3, and 
t/p*» corresponding to the good agreement of the asymp- 
tote to the second approximation with the exact t/p* 

curve . Since V^uT enters as a factor in all these equa- 
tions, there is obtained from the buckling determinant (38) 
independently of the value of u>, the following values 
for the coefficients in the expression for w: 



0 a = -0.380 0 X C 3 * -0.0470 0 % G 4 « -0.00565 0 1 
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From the comparison with the accurately computed coeffi- 
cients for u) == 172.5, it is seen that the buckling 
shape changes only slightly in the y direction for larger 
values of u) and is essentially distort ed only in the 
x direction, corresponding to equation (47) for the re- 
ciprocal wave length, 

The ratios in equation (31, x ) for U) = 172. 5 "become 
at larger values of a> and same central angles still 
smaller on account of the small change in the coefficients 
for the w expression with decreasing s s /b 2 . 

VI. SUMMARY ■ 



In the present paper the buckling condition for an 
infinitely long plate strip with constant curvature and 
simply supported longitudinal edges is set up for com- 
bined shear and longitudinal stress. Starting from the 
complete fundamental equations for the displacements a 
simple computation process is developed for plate strips 
of small curvature (central angle ^45°) with the aid of 
a Eitz expression. Por tlae case of pure shear stress the 

buckling condition is evaluated and for £ /.£ > 4.3 (u) > 

s v r 

40) a simple approximation formula derived. (See fig. 6 
and equation (49).) The computation results for the gen- 
eral case of combined shear and axial stresses will be 
reported on later. 
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Figs* 1.8,3,4,5,6,7 




Figure 2.- Displacement 

of middle 
surface of shell in the 
yz plane. 





n 


.A 




/ 


8 


3 






6 


f 


8 


3 




5 


6 


tef.f) 


f 




















*— 






8 




— 










— — 




mm 




mm 




3 






— 














— 
































S 


























e 




























f 


























2 


























3 


























¥■ 


























5 


























6 



























Figure 1.- Plate strip 
under 



Figure 5.- Scheme of coefficients of unknowns 
C13, Cg n in equation (81). 




Figure 4. -Shear 

forces 
on deformed plate 
element. 



Figure 6.- Critical 

shear 
stress of a curved 
plate strip as a 
function of 

\" 3 1— V* b*J 
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First approximation 
















— .Asymptote to second approximation 
For v » 0 at y ; 0 and y s h 










(third approxim 
to Leggett) . 
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For IT - 0 at y a 0 
and y a h. 
For v .a 0 at y « 0 
and y • b( according' 
to Leggett). 



Figure 3.« Section forces on curved 
plate element. 



Figure 7.- The reciprocal wave length 
£ for the smallest critical 
shear stresses as a function of v**/ 



